The long-wavelength behaviour of the exchange-correlation kernel in the Kohn-Sham 

theory of periodic systems 
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The polarization-dependence of the exchange-correlation (XC) energy functional of periodic in- 
sulators within Kohn-Sham (KS) density-functional theory requires a 0(l/q 2 ) divergence in the XC 
kernel for small vectors q. This behaviour, exemplified for a one-dimensional model semiconductor, 
is also observed when an insulator happens to be described as a KS metal, or vice- versa. Although it 
can occur in the exchange-only kernel, it is not found in the usual local, semi-local or even non-local 
approximations to KS theory. We also show that the test-charge and electronic definitions of the 
macroscopic dielectric constant differ from one another in exact KS theory, but are equivalent in 
the above-mentioned approximations. 
PACS numbers: 71.10.+X, 77.22. Ej 



I. INTRODUCTION 



We have recently reexamined |1|-|3|] the exchange- 
correlation potential, V^ c (r), that represents electronic 
exchange and correlation in the Kohn-Sham formulation 
of density- functional theory (DFT), in the case of macro- 
scopic bodies. Our analysis focused on periodic solids 
with a uniform macroscopic polarization: a crystal sub- 
jected to a homogeneous external electric field, or a polar 
crystal. We showed that the exact V xc inside the crystal 
is required to have an ultra-non-local dependence on the 
surface electron density, or, equivalently, on the macro- 
scopic polarization. In the present paper, we will show 
that the exchange-correlation kernel, K xc , is a power- 
ful vehicle for analyzing the presence of this dependence 
in the common approximations for exchange and corre- 
lation. The required wavevector-dependence of K xc is 
illustrated in the exact DFT of a model system. Links 
with the polarizability in insulators and metals will also 
be established. 

We first define some useful quantities and notation. 
Within the Kohn-Sham (KS) formulation of Density 
Functional Theory j| , a system of interacting electrons in 
its ground state, placed in an external potential, V ex t(r), 
is mapped onto a fictitious system of independent parti- 
cles in an effective potential 



^Ks(r) = Kxt(r) + Vk(r) + V xc (v), 



(1) 



under the requirement that both generate the same den- 
sity n(r). The combined external Kxt( r ) and Hartree 
Vfi(r) potentials give the electrostatic potential Vrc( r ) 
felt by a classical test charge (TC), while the exchange- 
correlation (XC) term, V KC (r), that subsumes all the ad- 
ditional many-body effects, acts only on KS electrons. 

The XC potential at point r is the first derivative of 
the XC energy with respect to the density at this point, 

Kc(r) = (2) 



The second derivative of this XC energy with respect to 
the densityis called the XC kernel: 



#xc(r,r') 



5n(r)Sn(r') 



(3) 



The explicit forms of E xc , V xc or K xc as functionals of the 
density are unfortunately unknown and calculations are 
usually performed within the local-density approxima- 
tion (LDA) , within semi-local approximations such as the 
generalized gradient approximations (GGA) or even 
within non-local approximations such as the weighted- 
density approximation (WDA) ||. 

In the same manner as V xc constitutes a sizeable con- 
tribution to the KS effective potential, K xc plays an im- 
portant role in the investigation of the responses of the 
KS system to static external pertubations. The density 
response of the interacting electron system to an external 
perturbation is described by the polarizability matrix x : 



5n(r)= / X {r,r')5V ext (r')dr'. 



(4) 



The so-called "proper part" of the density response, it [Q, 
relates the change of density to the change of the test- 
charge potential: 



Sn(r) = J 7r(r,r') <5P TC (r')dr'. 



(5) 



From the relation between TC and external potentials, 
these matrices are related by the following expression 
(making use of matrix notations): 



X 



V C , 



(6) 



where Vc stays for the Coulomb interaction. For periodic 
solids, the various quantities are most conveniently rep- 
resented in reciprocal space, where x(q+ G, q+ G'), for 
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example, may be thought of as a matrix in the reciprocal- 
space lattice vectors G and G', for each wavevector q, 
with the element G = G' = known as the head, the 
other elements in that row and column as the wings, and 
the rest of the matrix as the body p|. 

For the KS system, an independent-particle polar- 
izability matrix Xo is usually introduced, relating the 
change of density to that of the effective KS potential: 



5n(r) 



Xo(r,r')6V KS (r')dr'. 



(7) 



Using Eqs. ([!]) and (J^), this independent-particle po- 
larizability matrix is easily linked with the interacting- 
particle polarizability and its proper part: 



X- 1 + V C + K x 

1 +K* r , 



(8) 



Following the Adler and Wiser sum-over-state tech- 
nique , Xo can be directly computed from the KS wave- 
functions. As shown by Eq. M), the knowledge of K^ c is 
crucial in deducing x or n from Xo- This key role was 
recently recalled by Dal Corso, Baroni and Resta ]Tc| ]. 

In connection with the KS band-gap problem it was 
shown by Godby and Needs H ], within certain approx- 
imations, that the ground state of a periodic insulator 
is sometimes described as a metal in KS theory. This 
is particularly striking in that the long-wavelength be- 
haviour of the polarizability matrix x for an insulator is 
qualitatively different from that of a metal, while x is a 
ground-state quantity and should be correctly obtained 
within DFT. Godby and Needs emphasized that a highly 
nonanalytic K xc could allow such a phenomenon. 

Recently, we proposed that the exact XC energy func- 
tional is polarization-dependent in the case of periodic 
insulating solids submitted to a homogeneous electric 
field JlJ . After exploring in some detail the consequences 
of this finding, Aulbur, Jonsson and Wilkins showed 
the presence of a sizeable contribution of the polarization- 
dependence of E xc to the linear and non-linear optical 
response of real materials. For polar solids, careful treat- 
ment of the polarization is mandatory, even if no homo- 
geneous electric field is present M. Resta H tried to 
address the origin of these effects: he proposed that the 
polarization dependence should arise from the Coulomb 
coupling between the electron, in the bulk, and part of its 
correlation hole delocalized at the surface. At variance, 
Martin and Ortiz |l4|] , pointed out that the shape of the 
XC hole may already depend on the surface charge, and 
we demonstrated || that a polarization-dependence is ex- 
pected at the purely exchange level (without correlation) . 
In another paper, Martin and Ortiz ]l5[ | placed the new 
density-polarization functional theory in the perspective 
of important works of the seventies, and presented an 
alternative formulation of it. 



In Section II of the present paper, we link the 
polarization-dependence of the XC energy to a 0(1/ q 2 ) 
divergence of K xc in the limit of q — > 0, briefly sketched 
in Ref. |l(J, and illustrate this behaviour in the case of 
a one-dimensional semiconductor. An analysis of related 
points was given in Ref. jl2j . This requirement of the ex- 
act KS theory is not fulfilled in the usual approximations 
such as the LDA, the GGA, and the WDA, as shown 
in Section III. In Section IV, we observe that a 0(1/ q 2 ) 
behaviour of the exact XC kernel is able to resolve the 
"true insulator - KS metal" paradox. A unified treat- 
ment of K xc for insulators and metals is given in Section 
V, together with an analysis of the exchange-only kernel. 
We will finally discuss (Section VI) the consequences of 
these facts for the different definitions of the macroscopic 
dielectric constant. 



II. DIVERGENCE OF K X c 

Working on periodic insulators first, we adopt the same 
perturbative approach (long- wave method) as in Ref. , 
and briefly recall one of its central result. The change 
of external potential SV cx t produced by an infinitesimal 
homogeneous electric field S£ ex t is (written in one dimen- 
sion for brevity): 



SV ext (r) = lim SV cxt (q) (e^ r 

q^Q 

= lim simgr) 

9^0 q 



1 



(9) 



In reponse to this perturbation, the system will develop 
a change of density 8n(r). Owing to local field effects 
(Umklapp processes), it may contain contributions at dif- 
ferent (q + G) vectors (where G belongs to the reciprocal 
lattice). Within linear response, the long- wave part of Sn 
takes the form: 



Sn(r) 



lim q 8V sm(qr) 

q^Q 



(10) 



where 5V is the change of polarization for q = [JlTj . 

Generally, the self-consistent screening potential will 
also contain long-wave and more rapidly varying terms. 
Its long-wave part will include an Hartree contribution, 
corresponding to the screening of the applied field due to 
the Coulomb interaction. In Ref. Q|, we demonstrated 
that the polarization dependence of _E XC will manifest it- 
self through a homogeneous XC electric field 5£ xc , so that 
for q — > 0: 



5V xc (q) = 



2iq 



(11) 



From Eqs. (||) and (||), the XC kernel and potential are 
related by 5V XC = K xc Sn. Isolating in this equation the 
long-wave terms from the other contributions, we obtain, 
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in a generalized matrix notation (in which G stands for 
all non-zero vectors of the reciprocal lattice): 



1 5£_ 

q 2 

SV 



( -a ^ 
xc V 8n q +G 



(12) 



In order for the change in exchange-correlation field <5£ xc 
to be finite when a finite change of polarization SV takes 
place, the head of the exact exchange-correlation kernel 
matrix, K ]CC (q,q), must exhibit a 0(l/q 2 ) divergence in 
the limit of q -> 

The previous results can be illustrated in the case of 
a simple model one-dimensional semiconductor, already 
used in Ref . M , and for which K xc can be computed ex- 
actly (within the model). In this model, the sum of the 
external and Hartree potential is taken to be: 



Vext(x) + Vr(x) = V cos(27ra/a), 



(13) 



where a is the unit cell length. Moreover, a simple non- 
local self-energy operator is present, with the aim of mim- 
icking the relevant many-body effects. It has the same 
non-local form and same parameters as in Refs. pq,pj: 



E(r,r',u>) = 



/(*) + /(*') 



g(\x ~ x'\) 



(14) 



where f(x) = F a [1 — cos(27ra;/a)] is a negative function 
that has the cell periodicity and g(y) is a normalized 
gaussian of width w. We construct an exact KS theory 
for this model system by determining the local potential 
VkSj which, when filled with non-interacting electrons, 
reproduces the density obtained when including the self- 
energy operator [HE]- 

The proper part 7r of the polarizability of the model 
system can be obtained using the Adler and Wiser || 
sum-over-state technique applied to the eigenfunctions 
and eigenvalues of the reference Hamiltonian (that in- 
cludes S, which is taken to be independent of changes in 
the external potential). 

From our KS wavefunctions, that reproduce the 
same density, it is also possible to compute the KS 
independent-particle polarizability, Xoi from the Adler 
and Wiser technique. The relationship between tt and 
Xo is given by Eq. (0), so that we obtain : 



K xe = x; 1 -ir- 1 . 



(15) 



In Fig. [I] we have plotted the diagonal part of the com- 
puted K xc . The calculation was performed on a 80- unit- 
cell, which guarantees a convergence better than 0.7%. 
We observe a divergence in the limit of q — » . The inset 
exhibits its expected 0(1/ q 2 ) character. 

We note that the divergence of K xc has been obtained 
without including long-range correlation effects in our 
model semiconductor [[l9| . Inclusion of such effects would 
simply modify the coefficient of the K xc (q, q) divergence. 



III. APPROXIMATE FUNCTIONALS 

In the previous Section, we have linked the 
polarization-dependence of the XC energy with an 
0(1/ q 2 ) divergence of the XC kernel, and given an ex- 
ample of this behaviour. The link with a third concept, 
the ultra- non-local sensitivity of the exchange-correlation 
functional |Q was emphasized in Refs. &0]: a change in 
surface charge may have an influence on the bulk XC po- 
tential, independently of the distance between the point 
in the bulk and the surface. Any approximate XC func- 
tional can now be analyzed in the light of these charac- 
teristics of the exact functional. Although they have the 
same physical origin, each of them provides a different 
point of view. 

As regards the ultra-non-locality requirement, the be- 
haviour of LDA and GGA is clear: the corresponding 
exchange-correlation potentials at any point do not de- 
pend on the density outside the immediate neighborhood. 
The following analysis of the small wavevector depen- 
dence of the exchange-correlation kernel provides a more 
refined picture of the violation of this requirement. 

To start with, we consider the XC kernel in the LDA, 



K. 



LDA/ 



SV 



LDA 



Sn 



6(r-r'). 



(16) 



The Fourier transform of this kernel, diagonal in real 
space, is such that .?C^P A (q, q) is independent of q ppf . 
The gradient-corrected XC energy has the form 



K 



GGA r 



e xc [n(r), Vn(r)]dr, 



with the corresponding potential : 



V 



GGA 



(r) 



dn 



de x 



dVr 



(17) 



(18) 



The relationship between the long wavelength GGA-XC 
potential and the long wavelength density is therefore 
governed by the following kernel (compare with Eqs.(15) 
and (22) of Ref. EfJ): 



K. 



GGA 



(q,q) 



i 
o 



dn 2 



dv 



-2iVg Q / 
Jn 



dr 



In (dn)(d{d a n)) 

Ef d 2 e KC 
,/ a9lS Jn (d(d a n))(d(d n)) 



dr 



(19) 



The first term of the right member of this equation is 
the Fourier transform of the term appearing in the LDA. 
Eq. ([l9]) makes clear that the gradient corrections pro- 
duces terms with positive powers of q, but not the re- 
quired 0(l/q 2 ) divergence. 
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Contrary to what is expected by Mazin and Cohen [pi] , 
the GGA has therefore no apparent ability to improve 
upon the LDA behaviour in this respect. Mazin and Co- 
hen were also expecting the weighted-density approxi- 
mation ( WDA) , in which the XC functional is truly non- 
local, to improve upon the LDA behaviour. We now an- 
alyze this case. 

The exchange-correlation energy is given exactly by the 
adiabatic connection formula, in terms of the coupling- 
constant averaged pair-correlation hole, G(r,r'; [n]): 



SG hom (r,r';n(r; [n])) 8G hom {r, r'; n(r)) Sn(r; [n])) 



n(r) ^ r ' r/; [" ]) n(r / )drdr / (20) 



This exact expression was the starting point of the analy- 
sis of the XC hole by Resta @. In the WDA [§, the ex- 
act pair-correlation function at each point r is replaced by 
the pair-correlation evaluated for a homogeneous electron 
gas with "weighted" density, G hom (r, r'; n(r; [re])), where 
n(r; [re]) is determined from the density in the whole 
space so as to satisfy the requirement that the exchange- 
correlation hole integrates to — 1. The corresponding 
WDA-XC potential is made of three terms 

2 J |ro-r'| 
1 f ^ G h ° m (r,r ;n(r; [re])) 

+ 2j n{r) ^ 

1 ff n(r)re(r') 5G hom (r, r'; re(r; [re])) , 



r - r' 



Sn(r ) 



drdr'. (21) 



If the WDA-XC energy functional were ultra-non-local, 
arbitrarily remote surface charges would influence the 
XC potential, so as to induce the possibility of a linear 
XC potential in the bulk. The first or second terms of 
Eq. (H^) would present such a behaviour if the homoge- 
neous gas pair-correlation function G hom were replaced 
by some constant (independent of the distance between 
r and r'): in this case they would be similar to a Hartree 
potential, that exhibits the desired dependence upon sur- 
face charge. A partly delocalized exchange-correlation 
hole would lead to the same kind of behaviour [|l3| . But 
the pair-correlation for a homogeneous gas, G hom , decays 
on average as the inverse of the fifth power of the distance 
between the points M . This decay will make these terms 
insensitive to the surface charge. 

In the third term of Eq. (^l|), the derivative of the pair- 
correlation function with respect to some change of den- 
sity appears. Ortiz and Martin have interpreted the cor- 
responding term in the exact expression as coming from 
the polarizability of the exchange-correlation hole JM) , 
sensitive to the TC or KS homogeneous electric field. 
But we will see that the WDA exchange-correlation hole 
is not as sensitive as the exact exchange-correlation hole. 
The derivative of the pair-correlation function can be ob- 
tained thanks to the chain rule, 



£n(ro) Sn(r) 
The above-mentioned decay of G hc 



Sn(r ) 



(22) 



"" !n affects both fac- 
tors of the right-hand side. The first factor will exhibit 
the same decay as the homogeneous gas pair-correlation 
function, while the dependence of re(r; [re]) upon <5n(ro), 
driven by the requirement that the exchange-correlation 
hole integrates to — 1, will also be spatially short-ranged: 
a finite change of surface charge density on the surface of 
the crystal, times the r~ 5 decay and integrated over the 
whole surface, yields zero contribution in the thermody- 
namical limit. As a consequence, all the terms in Eq. ( |2l| ) 
are insensitive to arbitrarily remote surface charges, so 
that no polarization-dependence is present in V^ DA . 

Thus, in opposition to the suggestion of Mazin and 
Cohen ]2l] , |22| ] , we find that the possible non- locality 
of approximate functionals is not sufficient to gener- 
ate a polarization-dependence: one needs an ultra-non- 
local dependence. This might be attained by consider- 
ing either a model pair-correlation function that does 
not integrate to —1 in the bulk ]l3|] , or a polarizable 
exchange-correlation hole jlj] , as in the exact exchange 
approach M. 



IV. THE METAL - INSULATOR PARADOX 

Godby and Needs |ll]] observed that the KS band 
structure of semiconductors might present the charac- 
teristics of a metallic state (the absence of a band gap). 
We now argue that it is precisely the same 0(1/ q 2 ) be- 
haviour of K xc (q,q) that allows this "true insulator-KS 
metal" paradox to be understood. 

In order to make the analysis as simple as possible, 
we first ignore local field corrections, and also impose 
cubic symmetry. In this case, only the head of the dif- 
ferent matrices appearing in Eq. (^) must be taken into 
account, which means that this equation reduces to a 
scalar equality. The small wavevector behaviours of \ 
and Xo are well-known in both the metallic and the in- 
sulating case [Q . The independent-particle polarizability 
of a KS metallic ground-state behaves like 



\im Xo(q, (?) =7) 

q^O 

while for a cubic KS insulator, one has 
\im Xo(q,q) = otq 2 , 



(23) 



(24) 



where 7 and a are some negative constants. The head 

of the polarizability matrix %(g, q) for a metal (interact- 

2 

ing electrons, not KS electrons) behaves exactly as — |- 
in the long-wavelength limit (which corresponds to com- 
plete screening of the Coulomb potential), while for cu- 

2 

bic insulators, it is — where (3 is a positive constant 
smaller than one, describing the incomplete screening. 
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Now, one imposes no divergence in Xo (Si^z) (metal- 
lic KS ground state), and non-cancelling divergences of 
X (q>Q) an d Vc (incomplete screening of the insulating 
system), so that Eq. (ph without local fields becomes 



1 47T 47T TJr . , 

- = —3-0 + T + «), 



(25) 



which proves that K xc (q,q) must have a 0(l/q 2 ) diver- 
gence. This result establish a connection between the 
"true insulator-KS metal" paradox and the polarization 
dependence of the XC energy for insulators: in both 
cases, the small-wavevector behaviour of the XC kernel 
is similar. 



V. UNIFIED TREATMENT OF K X c FOR 
INSULATORS AND METALS 



If the Kohn-Sham system is insulating, x 1 ex ~ 
hibit the appropriate divergences. Following Ref. [|, 
the exchange-hole polarizes in a homogeneous KS elec- 
tric field, so that one obtains a polarization-dependence 
of the exchange-only kernel. On the other hand, if the 
Kohn-Sham system is metallic to start with, Xo 1 nas n0 
divergence, and V x (q) for short wavevectors will vanish. 

Thus, the corrections needed to obtain an insulating- 
like polarizability while the Kohn-Sham system is metal- 
lic are entirely due to the correlation kernel. In the 
adiabatic coupling constant construction of the XC en- 
ergy B , the system at different value of the coupling con- 
stant A will undergo a metal-insulator phase transition : 
at A = 0, the characteristic system is the Kohn-Sham 
metallic system, while at A = 1, one deals with the truly 
interacting, insulating, system. 



We now consider the treatment of insulators and met- 
als in a unique framework. By the way, we will also 
generalize the demonstration contained in the proceed- 
ing section to the case where local fields are included (we 
treat not only the head of K xc , Xo 1 an d tt -1 , but also 
their wings and body). 

Following the analysis by Pick, Martin and Cohen |t]] 
of Xo and 7r, we find thatXo 1 : for the KS non- interacting 
system, and 7r _1 , for the interacting system, have the fol- 
lowing, similar, non-analytic behaviour. In the insulat- 
ing case, the head of these inverse matrices diverges like 
0(l/q 2 ), the wing elements diverge like 0(l/q), while the 
body elements, though non-analytic, are non-divergent. 
In the metallic case, no element shows a divergence for 
small wavevectors. The divergences in K xc = x^, 1 ~ T > 
following Eq. (|i~5|) , are easy to deduce from these results, 
if we suppose that there is no fortuitous cancellation of 
divergences when the difference is taken. 

For the (usual) "true metal-KS metal" case, since 
there is no divergence in Xo 1 or ^ 1 no element of K xc 
need diverge for small wavevectors in order to correct 
the long-wavelength response (although separate non- 
analytic behaviour at non-zero wavevectors may be re- 
quired to accommodate a disparity between the true and 
KS Fermi surfaces). In all the other conceivable cases 
(the usual "true insulator-KS insulator" case, the ex- 
otic "true insulator-KS metal" case, and the hypotheti- 
cal "true metal-KS insulator" case), the head of K xc will 
diverge like 0(l/q 2 ), the wing elements will diverge like 
0(l/q), and the body elements will not diverge. 

Although the exchange-only kernel is able to generate 
the right divergences in the "true insulator-KS insulator" 
case ||, it is unable to do so in the "true insulator-KS 
metal" case. Indeed, the argument developed in Ref. ||] 
relies on the following result : 



^(q) = E^ 1 (q-q + G) 



5E X 



VI. THE MACROSCOPIC DIELECTRIC 
CONSTANT 



Now that the correct behaviour of K xc has been dis- 
cussed, we would like to investigate its consequences in 
the calculation of the macroscopic dielectric constant of 
insulators. At the macroscopic level, this quantity can 
be obtained as [E3f: 



1 + 4-7T 



dP 
88 



(27) 



where V is the macroscopic polarization and £ is the 
macroscopic electric field. 

In order to include the local-field effects in the compu- 
tation of Eryo , Adler and Wiser || relied on the dielectric 
matrix, relating effective and external potentials as fol- 
lows 

J2 e(q+G,q+G')£Kff(q+G')=<514xt(q + G), 

q+G' 

(28) 



and connected the macroscopic dielectric constant to the 
head of the inverse dielectric matrix: 

£00 = Mm — jr- -. (29) 



9^0 £ x (q>q) 



G 



<JVks(q+G) 



The demonstration of Eq. (|29|) was reported at the RPA 
level. The next step was to include correctly the correc- 
tion induced by the exchange-correlation effects . 

In the test-charge formulation of e pi[ , the effective 
potential appearing in Eq. (|2^) is chosen as the one ex- 
perienced by a hypothetical classical charge, <5Vtc = 
(26) bVc-xt, + SVu- Imposing the variations of the external po- 
tential <5K x t(q + G) for G 7^ to be zero, one deduces 
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= 1 - 



W H (q) 



£TC(^q) SV ext (q) + 5V H (q) 



= 1 - 



'q 2 SV cxt (q) + 5Vk(q) 



(30) 
(31) 



The relationship that exists in the long-wave approach 
between field and potential and between charge and po- 



larization allows one to recover Eq. (27) 



In the electron formulation, the effective potential ap- 
pearing in Eq. J2S|) is replaced by the one felt by the 
Kohn-Sham electrons, 5Vks — SV cxt + §Vn + SV XC . The 
dielectric constant is then : 



_L gVk(q) + ^xc(q) 

1 (q, q) SV ext (q) + S U H (q) + 6V XC (q 



(32) 



In the absence of polarization-dependence of E xc 
(SV xc (q) = 0), Eq. @ reduces to Eq. ©. Within 
the LDA, GGA or WDA, the test-charge and electron 
definitions of £oo are therefore identical. 

Going beyond these approximations, we observe that 
the electron definition, Eq. (|32]), does not reduce to 
the test-charge one, Eq. j30|), by a simple addition of 
the exchange-correlation contribution in the denominator 
(macroscopic field), but that the numerator must also in- 
clude a modified Coulomb interaction. In addition, it can 
be checked that the test-charge definition of the macro- 
scopic dielectric constant can be formulated as a second 
derivative of the (electric) free energy |25| with respect 
to a homogeneous electric field, which is not the case for 
the electron formulation. 

It is also possible to analyze the direct effect of the 
divergence of head and wing elements of K xc (q, q) on di- 
electric matrices. Following Singhal and Callaway [53 , 
we use the equality 



Sn = Xo[SV cxt + SVh + 6V XI 



(33) 



to find the form of the dielectric matrices in terms of %„, 
Vc and K xc : 



etc = 1 - VcXo [1 - K xc Xo] \ 



e e = 1 - VcXo - KxcXo- 



(34) 



(35) 



These expressions are governed by the products VcXo 
and K xc Xo- The head and wing elements of Xo behave 
like 0{q 2 ) and O(q) in the limit of q -> @|c]]. V c is 
diagonal, with a divergent 0(1/ q 2 ) head. In the small 
wavevector limit, their product, an asymmetric matrix, 
will have finite head and body, while the upper wing will 
diverge like 0(1/ q) and the lower wing will behave like 
0(q). If we completely ignore the effect of the exchange- 
correlation terms in Eqs. (|34|) and (p5|), the latter form 
is also the form of the dielectric matrices. 

If we take into account the effect of a non- divergent 
K xc , of the LDA type, the product -KxcXo will behave like 



Xo- The handling of such contributions in both Eqs. (|3J) 
and ( |35| ) will only affect the body of the dielectric ma- 
trices, in the small wavevector limit. By contrast, if we 
take into account the possible divergences of K xc , the 
product K xc Xo will behave like VcXo- As such, it will 
be able to affect the leading behaviour of all elements in 
both formulations of the dielectric matrix. 

The macroscopic dielectric constant being obtained 
from the head of the inverse dielectric matrix Eq. (|2S|), 
the presence of K xc will always affect its amplitude. How- 
ever, if one deals with a LDA- type K xc , its influence will 
be mediated by local fields only , while the true K xc will 
modify it more directly, in particular through a modifi- 
cation of the head of the dielectric matrix. 



VII. CONCLUSIONS 

In conclusion, we have seen that the polarization de- 
pendence of E xc imposes a condition on the form of the 
exchange-correlation kernel: its head must diverge like 
0(l/q 2 ) in the limit of q — > 0. This condition is not 
satisfied within the LDA, GGA and WDA. The absence 
of this divergence in the approximate functionals will af- 
fect the amplitude of the computed macroscopic dielec- 
tric constant by a finite amount. A similar divergence 
is also needed, in order to reproduce the correct dielec- 
tric response, when insulators are described as metals in 
Kohn-Sham theory. 
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